Existence of ℋ-matrix approximants to the inverse FE-matrix of elliptic operators with L∞-coefficients

نویسندگان

  • Mario Bebendorf
  • Wolfgang Hackbusch
چکیده

This article deals with the existence of blockwise low-rank approximants — so-called H-matrices — to inverses of FEM matrices in the case of uniformly elliptic operators with L∞-coefficients. Unlike operators arising from boundary element methods for which the H-matrix theory has been extensively developed, the inverses of these operators do not benefit from the smoothness of the kernel function. However, it will be shown that the corresponding Green functions can be approximated by degenerate functions giving rise to the existence of blockwise low-rank approximants of FEM inverses. Numerical examples confirm the correctness of our estimates. As a side-product we analyse the H-matrix property of the inverse of the FE mass matrix.

برای دانلود متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید

ثبت نام

اگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید

منابع مشابه

The spectral properties of differential operators with matrix coefficients on elliptic systems with boundary conditions

Let $$(Lv)(t)=sum^{n} _{i,j=1} (-1)^{j} d_{j} left( s^{2alpha}(t) b_{ij}(t) mu(t) d_{i}v(t)right),$$ be a non-selfadjoint differential operator on the Hilbert space $L_{2}(Omega)$ with Dirichlet-type boundary conditions. In continuing of papers [10-12], let the conditions made on the operator $ L$ be sufficiently more general than [11] and [12] as defined in Section $1$. In this paper, we estim...

متن کامل

Approximation of solution operators of elliptic partial differential equations by ℋ- and ℋ2-matrices

We investigate the problem of computing the inverses of stiffness matrices resulting from the finite element discretization of elliptic partial differential equations. Since the solution operators are non-local, the inverse matrices will in general be dense, therefore they cannot be represented by standard techniques. In this paper, we prove that these matrices can be approximated by Hand H2-ma...

متن کامل

Efficient inversion of the Galerkin matrix of general second-order elliptic operators with nonsmooth coefficients

This article deals with the efficient (approximate) inversion of finite element stiffness matrices of general second-order elliptic operators with L∞-coefficients. It will be shown that the inverse stiffness matrix can be approximated by hierarchical matrices (H-matrices). Furthermore, numerical results will demonstrate that it is possible to compute an approximate inverse with almost linear co...

متن کامل

ℋ-Matrix approximation for the operator exponential with applications

We develop a data-sparse and accurate approximation to parabolic solution operators in the case of a rather general elliptic part given by a strongly P-positive operator [4]. In the preceding papers [12]–[17], a class of matrices (H-matrices) has been analysed which are data-sparse and allow an approximate matrix arithmetic with almost linear complexity. In particular, the matrix-vector/matrixm...

متن کامل

Existence of ℋ-matrix approximants to the inverses of BEM matrices: The simple-layer operator

We consider the question of approximating the inverse W = V of the Galerkin stiffness matrix V obtained by discretizing the simple-layer operator V with piecewise constant functions. The block partitioning of W is assumed to satisfy any of the standard admissibility criteria that are employed in connection with clustering algorithms to approximate the discrete BEM operator V. We show that W can...

متن کامل

ذخیره در منابع من


  با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید

برای دانلود متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید

ثبت نام

اگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید

عنوان ژورنال:
  • Numerische Mathematik

دوره 95  شماره 

صفحات  -

تاریخ انتشار 2003